Physical Origin of the Universal Three-body Parameter in Atomic Efimov Physics 
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We address the microscopic origin of the universal three-body parameter that fixes the spectrum 
of few-atom systems in the Efimov regime. We identify it with a nonadiabatic deformation of the 
three-atom system, which explains its universal relationship with the van der Waals length scale 
observed in several experiments and numerical calculations with atoms, regardless of the depth of 
their interaction potential. 



In recent years, the investigation of Efimov physics [1 
12] , the universal physics of few particles interacting with 
nearly resonant short-range interactions, has developed 
tremendously, both on the experimental [3-22] and the- 
oretical fronts |23ff2"7] . The essence of this physics is 
the appearance of a universal 1/i? 2 attraction between 
three particles with an average separation R. This long- 
range three-body attraction, discovered by V. Efimov [T], 
emerges from the pairwise interactions, despite their fi- 
nite range b. It can be interpreted as an interaction be- 
tween two particles mediated by the third particle. Its 
strength is universally determined by the masses and 
quantum statistics of the particles. The Efimov attrac- 
tion extends from distances on the order of or larger than 
the range of the interaction b to distances on the order 
of \a\, where a is the scattering length of the pairwise in- 
teraction. It therefore requires \a\ > b, which is the case 
of resonant pairwise interactions. At the unitarity limit 
a — > oo, the Efimov attraction extends to infinite dis- 
tances. Decaying as 1/i? 2 , it supports an infinite number 
of bound states known as the Efimov trimers. Further- 
more, each state is related to the next excited state by a 
scale transformation, due to the scale invariance of 1/ R 2 
potentials [TJ [51] , so that the energy spectrum forms a 
geometric series. This constitutes the most remarkable 
and characteristic feature of the Efimov trimers. How- 
ever, since the three particles are attracted to each other, 
the physics at short separations comparable to the range 
b fixes the wave functions and spectrum of the Efimov 
trimers. 

Until recently, little has been known about this short- 
distance physics. Efimov's original investigation made 
use of the asymptotic two-body behaviour (or zero-range 
potential limit) to derive the three-body attraction, but 
could not address the short-distance region directly. Its 
effect on the longer-distance region was accounted for 
by a three-body boundary condition, expressed either as 
a phase in the three-body wavefunction or a log-periodic 
inverse length A known as the Efimov three-body param- 
eter pQ. In this long-distance picture, A may be regarded 
as the parameter that rcnormalises the low-energy effec- 
tive field theory [25]. This is the only known physical ex- 
ample of the renormalisation-group limit cycle [28 . Since 
the short-distance region involves the short-range details 
of the interaction potentials, A was long thought to be a 
non-universal quantity that is strongly dependent on the 



system. 

Later, it was found that A is universally determined 
in cases where a length scale larger than b arises in the 
problem, namely the case of a narrow Feshbach resonance 
in the pairwise interaction |29] , which entails a large and 
negative effective range setting the value of A |23j . and 
the case of particles with additional dipolar interactions, 
whose strength also sets the value of A (27]. In the ab- 
sence of such large length scales, however, it was believed 
that A ~ 1/b, but its precise value would vary by a factor 
spanning its log-period w 22.7 from one system to 

another, or even from one resonance to another within 
the same system |30| . 

However, several recent experiments with identical 
ultra-cold atoms [5H22] have revealed Efimov trimers 
and thereby determined their three-body parameters. 
In these experiments, rather broad Feshbach resonances 
are used, implying that the range b of the interactions 
between atoms is typically the van der Waals length 
ivdw = (mCs/H 2 ) 1 / 4 associated with the —Ce/r 6 tail 
of the open-channel potential |29j . The measured value 
of the three-body parameter expressed in units of £ vdw 
turned out to be surprisingly consistent between differ- 
ent atomic species [20], nuclear spin states [16], or even 
resonances of the same species |18| . 

A first attempt [3T| to explain these observations sug- 
gested that quantum reflection in the sum of pairwise 
— 1 /r 6 potentials could be the origin of such universality 
of the three-body parameter. This suggestion, however, 
was subsequently dismissed by a numerical study |32[ 
which uses model potentials to show that the universal- 
ity stems instead from an effective universal three-body 
repulsion that emerges at separations on the order of the 
van der Waals length. This three-body repulsion occurs 
whenever the pairwise interaction potential features a 
deep well supporting many two-body bound states, which 
is the case for the atomic species used in the experiments. 
It was therefore suggested that the deep well, by strongly 
accelerating two of the particles, thereby reducing their 
probability to be found at short distance, would provide 
a large kinetic energy to the three particles that would 
manifest itself three-body repulsion; however the 
physical mechanism behind this universal repulsion has 
remained unclear. On the other hand, another calcu- 
lation [3 lj with a realistic helium pairwise potential also 
found a three-body parameter consistent with the univer- 
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sal experimental results, even though the well of this po- 
tential is shallow and supports at most one bound state. 

In this letter, we elucidate the physics behind the 
universality of the three-body parameter in homonuclear 
atomic systems, which explains the results of both works 
|32| and pH] , and gives an intuitive understanding of the 
universal three-body parameter. 

Our starting point is the low-energy Faddeev equa- 
tion |37| . The three-body wave function in hyperspheri- 
cal coordinates reads 
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is the hyperradius, and ft 



designates all the angles and hyperangles. This is an 
exact decomposition over a basis of hyperangular wave 
functions $„(Sl;i?) normalised to 1. 

In the low-energy Faddeev approximation for three 
bosons, l>„ is reduced to V*, , „ „ ^"S a * :R ' . Here, a* = 

arctan ^ rj ^ is the Delves hyperangle in the ith Jacobi 
coordinate system {fjj~, Tijk)- By doing so, we ignore the 
dependence of <\> n on the directions of fjk and Pi jk, i-e- 
higher angular-momentum partial waves. This excludes 
accidental resonances with higher partial waves |30[ 132] . 
but otherwise gives a good description of three particles 
at low energy. 

The Faddeev component cf> n is chosen to be the eigen- 
solution with the eigenvalue A„ of the Faddeev equation 
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where m is the particle mass, V is the pairwise interaction 
potential, a mia = \j — a\, and a max = § — |-| — a\. We 
assume that three-body interactions are negligible |40| . 

Each solution <p n defines a channel n, and one finds 
in general that the hyperradial functions /„ (i?) are solu- 
tions of the coupled equations, 
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Figure 1: Three-body potentials U n (R) for different pairwise 
interactions at unitarity: soft van der Waals potential (blue) 
with up to 9 two-body bound states, Lennard- Jones poten- 
tial (green) with up to 9 two-body bound states, helium 
potential (red) rescaled to reach unitarity with two-body 
bound state. Note that only the case of the soft van der 
Waals potential with no bound state is significantly different 
from the other cases; all other curves nearly coincide and act 
as a universal repulsion that produces a universal three-body 
parameter. The dashed curve shows the asymptotic Efimov 
attraction. 



that P nn — due to the normalisation of <!>„, we arrive 
at a simple Schrodinger equation 
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with the three-body potential 
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This potential is a sum of an adiabatic (first term) and 
nonadiabatic (second term) contributions. The Efimov 
attraction manifests itself in this framework as the ap- 
pearance in a particular channel n — of a nega- 



tive eigenvalue Xq(R) 



-Sq at large hyperradii, with 



s fa 1.00624. One can show that Q 00 {R) -> 0{l/R 3 ), 
so that the potential Uq(R) tends to the l/R 2 Efimov 
attraction at large R. However, it is expected to grad- 
ually deviate from it at shorter distances comparable to 
the range of the interaction. 

We solved these equations numerically for two kinds of 
pairwise potentials with a van der Waals tail, namely the 
soft van der Waals potential and Lennard- Jones poten- 
tial: 



with the nonadiabatic couplings 



If we restrict our consideration to a particular channel 
n (neglecting couplings to other channels), and noting 
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where the length a is adjusted to produce a shape res- 
onance (divergence of the scattering length, leading to 
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Figure 2: Two-body probability density distribution \ip\ 2 (nor- 
malised to asymptote to 1) as a function of interparticle dis- 
tance for different two-body potentials: soft van der Waals 
potential (dotted blue) with up to 7 bound states, Lennard- 
Jones potential (solid green) with up to 7 bound states, 
helium potential (dashed red) rescaled to unitarity with 
bound state. The corresponding potentials are shown in faded 
colours (in units of the van der Waals energy ft 2 /mL 2 dw ). 



the unitarity limit a — > oo). There are several choices 
of c, corresponding to different potential well depths, or 
equivalently different numbers of two-body bound states. 
The resulting three-body potentials U n (R) for different 
potential depths are represented in Fig. [T] 

As in Ref. [32] , we find that in the limit of deep two- 
body potentials the three-body potential converges to a 
universal potential exhibiting a repulsion at short dis- 
tances. This universal short-range repulsion gives a uni- 
versal value to the three-body parameter of the Efimov 
trimers supported by this potential. For the Lennard- 
Jones potential, the three-body potential nearly coin- 
cides with the universal potential regardless of the num- 
ber of two-body bound states. Interestingly, the same is 
obtained with a realistic helium potential [33 rescaled 
to reach unitarity (34], which is qualitatively similar to 
a Lennard-Jones potential at unitarity with no bound 
state. This indicates that the apparent agreement be- 
tween the helium calculation in Ref. [34J and the exper- 
imental results |20_ is not a coincidence, but the conse- 
quence of a general universality class. 

To prove this, we need to understand the origin of 
this universal three-body repulsion. First, we argue that 
the universality originates from the fact that most of the 
physics happens in the van der Waals tail of the two-body 
potential. Such universality has been known for two par- 
ticles interacting with van der Waals interactions |351I36| . 
It is illustrated in Fig. [2] where the two-body probabil- 
ity density \(f\ 2 is represented for the different potentials. 
Within the tail of the potential, all densities nearly coin- 
cide. We now point out that even for a shallow potential 
with a short van der Waals tail, such as that of helium, 
the two-body density can be very similar to this univer- 
sal two-body density. The fact that it does not repro- 
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Figure 3: Three-dimensional contour plots of the probabil- 
ity distribution in Eq. ( 10 1 of finding a particle for a given 



separation of the two other particles (indicated by small gray 
balls connected by a black line). For clarity, the view is cut 
through the system along a plane containing the two particles. 
The top panels correspond to a separation of 3.0I/ vd w, while 
the bottom ones to a separation of 0.7L V dw- The darker, the 
denser the probability of finding the third particle. In order 
to appreciate the change of configuration between the panels, 
a typical location of the third particle is indicated by a small 
green ball connected to the two other particles by green lines. 
The left panels were computed from the zero-range Efimov 
theory at unitarity; they show the invariance of the Efimov 
configuration distribution with the size of the system. The 
right panels were computed for a Lennard-Jones pairwise po- 
tential at unitarity supporting three two-body bound states. 
At large separations, the probability is consistent with the 
Efimov configuration probability, but around each of the two 
particles there is a noticeable probability suppression region 
of spherical shape with a radius ~ i vdw . This suppression 
leads to an abrupt change of configuration probability when 
the particles come close. 



duce the short-distance oscillations does not bring any 
major difference, since the density is very small in this 
region. This suppression of density is due to the strong 
acceleration in the case of a deep soft-core van der Waals 
potential, or due to the hard-core repulsion in the case of 
helium or a shallow Lennard-Jones potential. Although 
of different physical origins, both cases lead to a small 
two-body density at short separations, so that most of 
the two-body physics happens in the universal van der 
Waals region. 

Next, we expound why this universality in the two- 
body sector leads to a repulsion for three particles. Our 
calculation shows that the repulsion originates from the 
nonadiabatic term Qoo m Eq. ([7]). This is a positive 
definite quantity, 
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This results in a significant gain of extra kinetic energy 
Qoo, creating the three-body repulsion. 

To check this point, we consider a simple model. We 
construct a trial hyperangular wavefunction which is the 
hyperangular function $o m the zero-range (Efimov) 
limit Eq. (Ill, multiplied by a Jastrow product of the 



universal two-body correlation ip causing the suppression 
of probability in the two-body sector, 



Figure 4: Comparison between the Faddeev three-body cal- 
culations (left) and the simple two-body correlation model 
described in the main text (right). The dashed curves show 
the nonadiabatic kinetic energy Qoo for Lennard- Jones poten- 
tials of different depths, corresponding to the unitarity limit 
with different numbers of two-body bound states ranging from 
to 4. The solid curves show the full three-body potential 
Uo(R) obtained by adding the adiabatic contribution Xo/R 2 
obtained from Faddeev calculations, which is shown by the 
dotted curves. 



owing to the normalisation of $o- This nonadiabatic ki- 
netic energy comes from a change of the hyperangular 
wavefunction $o with respect to the hyperradius. To vi- 
sualise this, we use $ to plot in Fig. |3] the probability 
density of finding a particle 3, for a given separation T\i 
of the two other particles 1 and 2, 



P{ri2, 3 ) = (sin2a 3 ) 2 
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In the zero-range limit (corresponding to the Efimov the- 
ory at large distance), <f> in Eq. (10) is given by 



00 (a) = sinh(s (- - a)) 



(11) 



which is independent of the hyperradius [ij. The prob- 
ability density therefore remains the same up to a scale 
transformation. In other words, the geometric config- 
uration of the three-particle system remains the same, 
namely the third particle is located close to one of the 
other two. The invariance of this Efimov configuration 
probability results in Qoo — 0- At shorter separations, 
however, this picture becomes inadequate because the in- 
teraction region is no longer neglig ible. In Fig. [3] one can 
clearly see two regions of suppressed probability near the 
two particles, for the reasons explained above. As the two 
particles come close, this suppression pushes out the Efi- 
mov configuration probability to a region forming a ring 
in-between the two particles. We argue that this change 
of configuration happens very suddenly, making it diffi- 
cult for the system to follow the Efimov channel adiabat- 
ically (see the animations in the supplemental materials) . 
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This simple ansatz leads to a very similar probability 
density. In particular, we have checked that it also leads 
to a sudden buildup of probability in the ring-shaped re- 
gion when two particles are close. One can also calculate 
the nonadiabatic kinetic energy Qoo from Eq. (pj. As 
expected, we find a sudden increase of this Qoo when the 
hyperradius is about £ vdw . As shown in Fig. ^] this re- 
produces very well the repulsion found using the solution 
of the Faddeev equation. This confirms our claim that 
the nonadiabatic change of configuration originates from 
an interplay of the suppression of two-body probability 
and the Efimov configuration. 

In conclusion, we have shown that the universality of 
the three-body parameter revealed in recent experiments 
with atoms and calculations is due to a three-body re- 
pulsion that originates from a nonadiabatic deformation 
of the three-atom system. This deformation is caused 
by the suppression of two-body probability at short dis- 
tance, which turns out to be universal for atomic systems, 
even for shallow potential species such as helium. This 
makes the three-body parameter universal in these sys- 
tems. More generally, we expect to find such universality 
in any class of systems where the two-body suppression 
has a universal behaviour. 

Finally, we point out that this universality is likely to 
extend to more than three particles. As the nonadia- 
batic deformation prevents three particles to come close, 
it is natural to expect that more than three particles 
cannot come close either. Therefore, we conjecture that 
iV-body atomic systems at low energy are universally de- 
termined by the two-body sector information only: the 
s-wave scattering length and the van der Waals length. 
This could explain the numerical evidence of universality 
obtained in previous works using realistic atomic inter- 
actions [381 I3U] . It would also remove the necessity of 
a four-body parameter in these systems, which has been 
advocated on different grounds in the scaling limit of ex- 
cited Efimov states 1251 1251. 
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This is justified to the extent that the range of the phys- 
ical three-body force does not exceed the range of the 
induced three-body repulsion found in this study. 



